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I Abstract. Three members of the new class of pulsating sdB 

stars (sdBV or EC 14026 stars) are analysed from Keck HIRES 
spectra using line blanketed NLTE and LTE model atmo- 
spheres. Atmospheric parameters (Tcff, log g, log(He/H)), 
metal abundances and rotational velocities are determined. A 
careful investigation of several temperature indicators, i.e. line 
profile fitting of Balmer and helium lines, the ionization equi- 
libria of helium, nitrogen and silicon gave consistent results for 
Feige48 and KPD2109+4401 to within a few hundred Kelvin. 
However, for PG 1219H-534 considerably higher effective tem- 
perature estimates were derived from the ionization equilib- 
ria of nitrogen (36 800K) and helium (34 400 K) than from 
the Balmer line profile fitting (33 200 K). A systematic differ- 
ence in the gravity derived from NLTE and LTE models was 
observed, the NLTE gravities being slightly lower, by up to 
0.1 dex, than the LTE results. 

As is typical for sdB stars, all programme stars are found 
to be helium deficient, with a helium abundance ranging from 
1/80 solar for Feige48 to 1/3 solar for PG1219H-534, prob- 
ably due to diffusion. Most metals are also depleted. The 
abundances of C, O, Ne, Mg, Al and Si in the high gravity 
programme stars KPD2109H-4401 and PG 1219H-534 are con- 
siderably lower than in the lower gravity stars Feige48 and 
PG 1605H-072 (Heber et al., 1999) which could be explained 
by an equilibrium between gravitational settling and radiative 
levitation. Surprisingly iron is solar to within error limits in all 
programme stars irrespective of their gravity, confirming pre- 
dictions from diffusion calculations of Charpinet et al. (1997). 

The metal lines are very sharp and allow the microtur- 
bulent velocity to be constrained to be lower than 5km/s 
(KPD 2109H-4401, PG 1219H-534). Also the projected rotational 
velocities have to be very low (v,ot sini < lOkm/s). For 
Feige48 the limits are even tighter {\ micro <3km/s, v^ot sini 
< 5 km/s). 



Send offprint requests to: U. Heber 

* Based on observations obtained at the W.M. Keck Observatory, 
which is operated by the Califomian Association for Research in As- 
tronomy for the California Institute of Technology and the University 
of California 



1. Introduction 

It is now well established that the hot subluminous B stars can 
be identified with models of the extreme Horizontal Branch 
(EHB, Heber, 1986, Saffer et al. (1994). Their photospheric 
chemical composition is governed by diffusion processes lead- 
ing to strong peculiarities (e.g. deficiencies of helium, carbon 
and silicon in some stars) with strong variations from star to 
star (for a review see Heber, 1998). 

Recently, several sdB stars have been found to be pulsating 
(termed EC14026 stars after the prototype, see O'Donoghue et 
al. 1999 for a review), defining a new instability strip in the HR- 
diagram. The study of these pulsators offers the possibility of 
applying the tools of asteroseismology to investigate the struc- 
ture of sdB stars. The existence of pulsating sdB stars was pre- 
dicted by Charpinet et al. (1996), who uncovered an efficient 
driving mechanism due to an opacity bump associated with 
iron ionization in EHB models. However, in order to drive the 
pulsations, iron needed to be enhanced in the appropriate sub- 
photospheric layers, possibly due to diffusion. Subsequently, 
Charpinet et al. (1997) confirmed this assumption by detailed 
diffusion calculations. Even more encouraging was the agree- 
ment of the observed and predicted instability strip. 

Eighteen pulsating sdB stars are well-studied photometri- 
cally (see O'Donoghue et al. 1999 for a review, and Billeres et 
al. 2000, Silvotti et al. 2000 and 0stensen et al. 2000 for more 
recent discoveries). A precise knowledge of effective tempera- 
ture, gravity, element abundances and rotation is a prerequisite 
for the asteroseismological investigation. 

We selected four EC 14026 stars for a detailed quantitative 
spectral analysis (see Fig. Il]): PG 1605+072 was chosen be- 
cause it has the lowest gravity and, therefore, has probably 
already evolved beyond the extreme horizontal branch phase. 
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Fig. 1. Position of the programme stars in the (Toff, logg) plane 
and comparison with other EC 14024 stars which are well- 
studied photometrically. The position of the zero age extreme 
Horizontal Branch (ZAEHB), the terminal age extreme Hori- 
zontal Branch (TAEHB), the helium main sequence and evo- 
lutionary tracks for EHB stars (Dorman et al., 1995) are also 
shown. 



It also displays the richest frequency spectrum amongst the 
EC 14026 stars (>50 periods have been identified, Kilkenny 
al. 1999). The results of an spectral analysis have already been 
reported by Heber et al. (1999). Recently, Kawaler (1999) pre- 
dicted from his modelling of the pulsations that PG 1605+072 
should be rotating. This prediction was confirmed by the analy- 
sis of optical spectra from which a projected rotational velocity 
of Vrot sini=39km/s was derived (Heber et al., 1999). In this 
paper we present the spectral analysis of three additional pul- 
sating sdB stars. PG 1219+534 was chosen because it has the 
shortest pulsation periods and has a helium abundance larger 
than most other sdB stars (O'Donoghue et al., 1999). Feige48 
was selected because it is the coolest of all EC 14026 stars 
known (Koen et al., 1998). The effective temperature and grav- 
ity of KPD2109H-4401 (Billeres et al., 1998) places the star in 
the parameter space where most EC 14026 stars are found (see 
Fig. 1). For Feige48 and KPD 2109+4401 only few frequen- 
cies (four for the former, five for the latter) have been found so 
far. 



2. Observations 

High resolution optical spectra were obtained with the HIRES 
echelle spectrograph (Vogt et al. 1994) on the Keck I telescope 
on July 20, 1998 using the blue cross disperser to cover the full 



wavelength region between 3700A and 5200A at a resolution 
0.09A. 

The spectra are integrated over one pulsation cycle or more 
since the exposure times (600-900s) were long compared to 
the pulsational periods. 

The standard data reduction as described by Zuckerman & 
Reid (1998) resulted in spectral orders that have a somewhat 
wavy continuum. In order to remove the waviness we used the 
spectrum of HI 504+65 (a very hot pre-white dwarf devoid of 
hydrogen and helium, Werner 1991) which was observed in the 
same night. Its spectrum has only few weak lines of highly ion- 
ized metals in the blue (3600^480A) where the strong Balmer 
lines are found in the sdB stars. Therefore we normalized indi- 
vidual spectral orders 1 to 20 (3600-4480A) of the sdB stars by 
dividing through the smoothed spectrum of HI 504+65. The re- 
maining orders were normalized by fitting the continuum with 
spline functions (interpolated for orders 26 and 27 which con- 
tain H/3). Judged from the match of line profiles in the over- 
lapping parts of neighboring orders this procedure worked ex- 
tremely well. Atmospheric parameters determined from indi- 
vidual Balmer lines are found to be consistent with each other 
except for H/3. Since this might be a result of interpolation er- 
rors in the normalization process, we excluded H/3 from the fit 
procedure. 

Besides the strong Balmer lines of hydrogen, helium lines 
are present in all three stars. However, He II 4686A is detected 
in KPD2109H-4401 and PG 1219H-534, only. In addition, weak 
photospheric metal lines can be identified in the spectra of all 
programme stars. However, the number of detectable photo- 
spheric lines differs considerably. The largest number of metal 
fines is present in Feige48 (C, N, O, Ne, Mg, Si, Al, S and 
Fe). In PG 1219+534 and KPD2109H-4401 only N, S and Fe 
are detectable. Also interstellar Ca II lines are present and are 
found to have a complex structure (three to four components, 
see appendix). 



3. Spectral analysis 

The simultaneous fitting of Balmer line profiles by a grid of 
synthetic spectra has become the standard technique to de- 
termine the atmospheric parameters of hot high gravity stars 
(Bergeron et al. 1992). The procedure has been extended to in- 
clude helium line profiles as well and applied successfully to 
sdB stars by Saffer et al. 1994). It will be referred to as Saf- 
fer's procedure throughout this paper The Balmer lines (H7 to 
H 12), He I (4471 A, 4026 A, 4922A, 47 13 A, 50 16 A, 5048 A) 
and He II 4686A lines are fitted to derive all three parameters 
simultaneously. 

The analysis is based on an updated version of the LTE 
model atmosphere code of Heber et al. (1984a) which includes 
metal line blanketing using Kurucz' ATLAS6 Opacity Dis- 
tribution Functions. A large grid of models is calculated for 
various helium abundances and solar metalicity as well as for 
metal poor ([M/H]=-2.0) composition (see Heber et al. 1999b). 
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Lemke's versionjj of the LINFOR program (developed origi- 
nally by Holweger, Steffen, and Steenbock at Kiel University) 
is used to compute a grid of theoretical spectra which include 
the Balmer lines Hq, to H22 and He I and He II lines. As can 
be seen from Tables nl and @ the metalicity dependence of 
the results is marginal, the effective temperature changes by 
600 K, the gravity by 0.02 dex and the helium abundance by 
0.09 dex at the most when we apply Saffer's procedure with 
the metal poor LTE grid ([M/H]=-2.0) instead of the solar com- 
position grid. Tests using models calculated with Kurucz' AT- 
LAS9 code were also performed. No significant deviations of 
the fit results from those listed in Table 1 were found. 

In order to investigate the role of NLTE effects we repeated 
the analyses using a grid of H-He line blanketed, metal free 
NLTE model atmospheres (Napiwotzki 1997), calculated with 
the ALI code of Werner & Dreizler (1999). Applying Saffer's 
procedure with the NLTE model grid (see Tables [I], ^ and g) 
yields T^ff and log ^^^^^ almost identical to that obtained from 
the LTE grid. We therefore conclude that the Balmer and he- 
lium lines are not vulnerable to NLTE effects. 

However, a systematic difference in log g persists, the LTE 
values being higher by 0. 1 dex than the NLTE results (see Ta- 
ble nl). Since its origin is obscure, we finally adopted the aver- 
aged atmospheric parameters given in Tables nl B^and pi 

The metal lines are sufficiently isolated to derive abun- 
dances from their equivalent widths again using the LINFOR 
program. Oscillator strengths were taken from the critical com- 
pilation of Wiese et al. (1996) for C, N and O, from the Opac- 
ity Project (Seaton, 1987) for Nell Mgll, and AlIII using 
the Opacity Project data base "TOPbase" (Cunto & Mendoza, 
1992) at the CDS (Strasbourg), from Becker & Butler (1990) 
for Si, Wiese et al. (1969) for S III, and Ekberg (1993) and Ku- 
rucz (1992) for Felll. The oscillator strengths for some lines 
of light elements which could not be found in the papers cited 
were supplemented from Kurucz' line list, from which also the 
damping constants for all metal lines were extracted. 

Equivalent widths of metal lines between 5mA and 79mA 
were measured. Since fewer metal lines are present in the 
spectra of KPD 2109H-4401 and PG 1219H-534 than in those of 
Feige48 and PG 1605H-072, we use the latter spectra as guid- 
ance to determine upper limits for metal abundances of the for- 
mer from the absence of the presumedly strongest lines of the 
ions in question. An upper limit to their equivalent width of 
5mA was assumed. 

3.1. Feige48 

Since no He II line can be detected, the helium ionization equi- 
Hbrium can not be evaluated. 

The results of the analysis of the H and He I lines for Feige 
48 are listed in Table n^ and compared to published values ob- 
tained from low resolution spectra. The formal errors of our 
fits are much smaller than the systematic errors (see below). 



Table 1. Atmospheric parameters for Feige 48 from different 
methods, see text 



method 


Tcff[K] logg logSJ^ 


Koen et al. (1998) 
LTE: H+He [M/H]=0.0 
LTE: H+He [M/H]=-2.0 
NLTE: H+He 


28 900±300 5.45±0.05 - 

29 500 5.53 -2.94 
30 100 5.53 -2.90 
29 600 5.43 -2.93 


adopted 


29 500±300 5.5 ±0.05 -2.93±0.05 



Table 2. Metal abundances for Feige 48 compared to solar 
composition, n is the number of spectral lines per ion. 



ion 


n 


lege 


[M/H] 


CHI 


2 


7.36±0.03 


-1.27 


NH 


16 


7.47±0.17 


-0.48 


NHI 


2 


7.80±0.11 


-0.15 


OH 


29 


7.79±0.12 


-1.09 


Nell 


3 


7.10±0.31 


-0.69 


MgH 


1 


6.91 


-0.6 


AlIH 


2 


5.50±0.18 


-0.89 


Si HI 


4 


6.39±0.32 


-1.01 


Si IV 


1 


6.17 


-1.07 


SHI 


1 


6.15 


-1.01 


FelH 


26 


7.55±0.19 


+0.13 



For a description see 
erlangen.de/~ai26/linfit/linfor.html 



http://a400.sternwarte.uni- 



The fit of the lines using the NLTE model grid is displayed 
in Fig. g. The agreement with the results from low resolution 
spectra analysed with similar LTE model atmospheres (Koen et 
al. 1998) is encouraging. 

Helium is deficient by a factor of 80 and therefore Feige 48 
has the lowest helium abundance among our programme stars. 

Three species are represented by two stages of ionization 
(C II and C III, N II and N III, Si III and Si IV). Since these 
line ratios are very temperature sensitive at the temperatures 
in question, we alternatively can derive Tes and abundances 
by matching these ionization equilibria. Gravity is derived sub- 
sequently from the Balmer lines by keeping Toff and log ^iiii 
fixed. These two steps are iterated until consistency is reached. 
C II is represented by the 4267A line only, which is known to 
give notoriously too low carbon abundances. Indeed the car- 
bon ionization equilibrium can not be matched at any reason- 
able Toff. The ionization equilibrium of N indicates a slightly 
higher effective temperature of rcff«30700K than derived 
from the Balmer and neutral helium lines using Saffer's pro- 
cedure whereas the silicon ionization equilibrium indicates a 
slightly lower effective temperature of Toff «28 000 K. How- 
ever, in the case of the nitrogen ionization equilibrium, N III 
is represented by two weak lines only, and in the case of sil- 
icon Si IV is represented by one weak line only. Considering 
the measurement errors for these lines, the nitrogen and silicon 
ionization equilibria can be considered to be in good agreement 
with the fit of the Balmer and neutral helium lines. 

The effective temperature derived from the metal ionization 
equilibria is higher than that from Saffer's procedure. Whether 
this may be caused by an NLTE effect has to be investigated 
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Fig. 2. B aimer and He line profile fits for Feige 48 of the HIRES spectrum. 



from detailed NLTE calculations for N and Si which are beyond 
the scope of this paper. 



Although several O lines are available in the spectrum of 
Feige 48, it was impossible to determine the microturbulent ve- 
locity in the usual way, i.e. by removing the slope of the O 
abundances with equivalent widths, due to the lack of suffi- 
ciently strong lines. However, the spectrum is very sharp lined 
and several close blends are resolved, e.g. the doublet struc- 
ture of Mg II 448 1 A is clearly seen (Fig. ||). The Mg II doublet 
is best fitted with no extra broadening (see Fig. Bl top panel), 
i.e. Vmicro <3km/s. Alternatively, we can constrain the pro- 
jected rotational velocity (see Fig. pi bottom panel), an upper 
limit of Vrot sin i<5 km/s was adopted. Microturbulent and pro- 
jected rotational velocities can also be constrained from un- 
blended lines such as Si III 4552, 4567A although the limits 
are not as strict {v micro <5km/s, Vrot sin i< 10 km/s) as those 
derived from the Mgll doublet. We adopt v„iicro=Okm/s. An 
increase to 5 km/s translates into small systematic abundance 
uncertainties of 0.05dex for most ions. A temperature uncer- 
tainty of ATofr=1000 K translates into abundance uncertainties 
of less than 0. 1 dex. Hence systematic errors are smaller for 
most ions than the statistical errors. The resulting abundances 
are Hsted in Table g and plotted in Fig. ^ 




4X/i 

Fig. 3. Fit of the Mg II doublet for microturbulent velocities 
(top) of Vmicro=0, 5 (dotted), and 10 km/s (dashed) and pro- 
jected rotational velocities (bottom) of Vrot sini=0, 5 (dotted), 
and 7 km/s (dashed) 

3.2. KPD 2109+4401 

As in the case of PG 1605+072 (paper I), Hell 4686A can 
be measured in the spectrum of KPD 2109+4401 allowing 
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Table 3. Atmospheric parameters for KPD 2109+4401 from 
different methods, see text 



Table 5. Atmospheric parameters for PG 1219+534 from dif- 
ferent methods, see text 



method 


Tcft [K] 


logg 


log "H- 


Billeres et al. (1998) 


31200 


5.84 




LTE: H+He [M/H]=0.0 


31200 


5.77 


-2.22 


LTE: H+He [M/H]=-2.0 


31400 


5.79 


-2.31 


LTE: He 


32400 


5.81 


-2.13 


NLTE: H+He 


31200 


5.69 


-2.30 


NLTE: He 


32 500 


5.73 


-2.18 


adopted 


31800±600 


5.76±0.05 


-2.23±0.1 



method 


Teff [K] 


logg 


log "Ho 

& riH 


Koen et al. (1999) 


32 800±300 


5.76±0.04 




LTE: H+He [M/H]=0.0 


33 200 


5.93 


-1.60 


LTE: H+He [M/H]=-2.0 


33 200 


5.95 


-1.62 


LTE: He 


35 200 


6.03 


-1.41 


NLTE: H+He 


33 300 


5.85 


-1.56 


NLTE: He 


35 400 


5.87 


-1.38 


adopted (see text) 


34 300+^""" 

J+JUU_ 1000 


5.95±0.1 


-1.5±0.1 



Table 4. Metal abundances for KPD 2109+4401 compared to 
solar composition, n is the number of spectral lines per ion. 



ion 


n 


loge 


[M/H] 


CIH 


- 


<6.3 


< -2.3 


Nil 


2 


7.26±0.30 


-0.69 


NIH 


2 


7.23±0.03 


-0.72 


on 


- 


<6.6 


< -2.3 


NeH 


- 


<6.5 


< -1.2 


Mgll 


- 


<5.9 


< -1.6 


AlHI 


- 


<5.3 


< -1.1 


Si III 


- 


<5.4 


< -2.1 


Si IV 


- 


<5.4 


< -2.1 


SHI 


2 


6.50±0.02 


-0.66 


FeHI 


12 


7.63±0.2I 


+0.21 



Table 6. Metal abundances for PG 12 19+534 compared to solar 
composition, n is the number of spectral lines per ion. 



ion 


n 


loge 


[M/H] 


CIH 


- 


<6.6 


< -2.0 


Nil 


7 


7.79±0.08 


-0.16 


NIH 


5 


7.79±0.16 


-0.16 


on 


- 


<7.0 


< -1.9 


NeH 


- 


<6.7 


< -1.0 


Mgll 


- 


<6.2 


< -1.3 


AlHI 


- 


<5.7 


< -0.7 


Si IV 


- 


<5.6 


< -1.9 


SHI 


2 


6.86±0.06 


-0.30 


FeHI 


2 


7.56±0.16 


+0.13 



the heUum ionization equilibrium to be exploited as well 
as Saffer's procedure. Both methods give parameters in res- 
onable agreement (see Table 0) and the averaged numbers 
(Teff=31 800K, log g=5.79 log(He/H)= -2.22) were adopted 
for KPD 2109+4401. The hne profiles calculated from this 
model reproduces the observed spectrum well (see Fig.Q). The 
derived Toff and log g are also in good agreement with the re- 
sults from low resolution spectroscopy (Billeres et al., 1998). 
Nitrogen is also present in two stages of ionization (N II and 
NIIl) and the effective temperature derived from its ioniza- 
tion equilibrium is in perfect agreement with the results from 
Balmer and helium line fitting. 

Helium is depleted by a factor of 17 with respect to the sun. 

From metal line profile fitting, as described for Feige 48, 
the projected rotational velocity of KPD 2109+4401 is con- 
strained by our spectra to v sini < 10 km/s. The micro turbulent 
velocity is constrained to Vmicro <5 km/s. Again we adopted 
Vmjcro=0km/s for further analysis. 

Abundances of nitrogen, sulfur and iron were derived from 
the equivalent widths of their spectral lines. Upper limits for the 
abundances of carbon, oxygen, neon, magnesium, aluminium, 
and silicon were derived from the absence of the presumedly 
strongest lines of those ions. 

The resulting abundances are listed in Table Q and plotted 
in Fig. H 



Fontaine (1999, priv. comm.) obtained an improved light 
curve of KPD 2109+4401 and detected some structure in the 
Fourier spectrum, in particular, a symmetric quintuplet, which 
may be due to rotational splitting. A rotational period of about 
6.6 days (or longer if not all the components of an n-tuplet are 
seen) has been inferred. Assuming a mass of O.5M0 we derive 
a radius of 0. 16i?0 from its gravity, which results in a rotational 
velocity of 1.2 km/s, too small to be detected in our HIRES 
spectra. 

3.3. PG 1219+534 

As already pointed out by Koen et al. (1999) the helium 
lines are stronger than in most other sdB stars. He II 4686A 
is detected allowing the He ionization equilibrium to be 
used as a temperature indicator Tcff and log g derived with 
Saffer's procedure agree well with results from low res- 
olution spectroscopy (Koen et al., 1999). However, unlike 
PG 1605+072 and KPD 2109+4401, the helium ionization 
equilibrium and Saffer's procedure yield discrepant results (see 
Table |): Teff=33 200K, log g=5.93, log(He/H)=-1.60 (Saffer's 
procedure. Fig. pi top panel) and Toff=35 200K, log g=6.03, 
log(He/H)=-1.41 (He ionization equilibrium. Fig. Bl bottom 
panel). At the lower Toff the Balmer lines are well matched 
throughout the entire profile, whereas for He II 4686A there 
is a significant mismatch (see Fig. q|, top panel). At the higher 
TofT He II 4686A is well reproduced, but the Balmer line cores 
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Fig. 4. Balmer and He line profile fits for KPD 2109+4401 of the HIRES spectrum. 



are not reproduced at all (see Fig. B, bottom panel). The line 
cores of He I 4026 A and 4471 A cannot be reproduced by ei- 
ther model. Despite of its high gravity PG 1219+534 has an 
unusually high helium abundance, i.e. helium is deficient by a 
factor of 2 to 5, only. 

Besides the Balmer and the helium lines, metal lines of 
N II, N III, S III and Fe III are present. The ionization equi- 
librium of nitrogen (7 N II and 5 N III lines) can be used to 
determine the effective temperature, which requires the effec- 
tive temperature to be as high as T'cff=36800K. At the low Tcs 
indicated by the Balmer line fit, the nitrogen abundance derived 
from N II lines differs from that derived from N III lines by al- 
most 1 dex. 

The abundances were calculated from the model most con- 
sistent with the nitrogen ionization equilibrium and results are 
listed in Table and plotted in Fig. pi The uncertainty in Tcs 
was taken into account when we estimated the upper limits for 
carbon, oxygen, neon, magnesium, aluminium and silicon from 
the absence of their presumedly strongest lines. 

Since the metal lines are very sharp we can constrain either 
the microturbulent velocity or the projected rotational velocity 
as demonstrated for Feige 48. Due to the weakness of the lines 
the limits are not as tight as for Feige 48: v micro <5 km/s, v sini 
< 10 km/s. 



4. Discussion 

We have determined atmospheric parameters (Tcs, log g, 
He/H), metal abundances and projected rotation velocities from 
time averaged, high resolution spectra for four pulsating sdB 
stars (Paper I, this paper). The effective temperatures and grav- 
ities (see Fig. 1) confirm that KPD2109H-4401, PG 1219H-534 
and Feige 48 are bona-fide extended horizontal branch stars, 
the latter being on the terminal EHB line, whereas the gravity 
of PG 1605H-072 is too low to be consistent with an EHB nature 
and might indicate a post- EHB evolutionary status for this star 

4.1. Atmospheric parameters 

Using line blanketed LTE model atmospheres for solar and 
metal poor composition ([M/H]=-2.0) as well as line blanketed 
zero metalicity NLTE model atmospheres we found that 

- metalicity effects are small. 

- NLTE effects are unimportant for the temperature determi- 
nation. 

- there is an offset between gravities determined from the 
NLTE grid and the LTE grid. The NLTE gravities are 
slightly lower, by up to 0. 1 dex. 

A careful investigation of several temperature indicators, 
i.e. line profile fitting using Saffer's procedure, the ionization 
equilibria of helium, nitrogen and silicon gave consistent re- 
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Fig. 5. Top: Balmer and He line profile fits for PG 1219+534 of the HIRES spectrum. Note the mismatch of the He II 4686 A line 
profile and the cores of He 1 4026 A and 4471 A. 

Bottom: He line profile fits for PG 1219+534 of the HIRES spectrum to determine Tcff and log(He/H) simultaneously, log g is 
adjusted to match the Balmer line wings. Note the mismatch of the cores of the Balmer lines and of He 1 4026A and 4471 A. 
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suits for Feige48, KPD2109+4401 and PG 1605+072 (paper 
I). However, for PG 1219+534 a considerably higher effective 
temperature was derived from the ionization equilibria of nitro- 
gen (36 800K) and hehum (34 400 K) than from the Hne pro- 
file fitting (Saffer's method, 33 200 K). What causes these dis- 
crepant results ? NLTE effects might be considered. However, 
our NLTE calculations show that this is not the case for helium. 
For nitrogen such calculations are beyond the scope of this pa- 
per. In view of the consistent results for the other stars, NLTE 
effects for the N ionization equilibrium in PG 1219+534 appear 
unlikely. 

PG 1219+534 has the highest He abundance and the short- 
est pulsation period of our programme stars. Hence the dis- 
crepancy could be related to these properties. Since the helium 
lines are stronger than in the other programme stars, line broad- 
ening is more important for the He I lines in the spectrum of 
PG 1219+534. We used tabulations of Barnard et al. (1969), 
Shamey (1969), Barnard et al. (1974, 1975) and Griem (1974) 
to compute the He I line profiles. New calculations have be- 
come available recently (Beauchamp et al., 1997). However, 
these tables are tailored for the use in white dwarf atmospheres 
and cannot be used for sdB stars because they do not extend 
to the significantly lower densities in the outer layers of sdB 
atmospheres. An extension of Beauchamp et al.'s tables would 
be useful for the sdB star line profile synthesis. For He II lines 
we used the tabulations of Schoning & Butler (1989). 

Due to diffusion the helium distribution in the atmospheres 
could be inhomogeneous. If helium accumulates in the deeper 
atmospheric layers, the predicted He I line profiles, however, 
would be even broader and shallower, increasing the discrep- 
ancy with observations. A vertical stratification of helium, 
therefore, appears to be unlikely. 

On the other hand the short pulsational periods 
(P=128. Is. ..148.8s) observed for PG 1219+534 might 
indicate that the model assumption of hydrostatic equilibrium 
does no longer hold for the outer layers of the atmosphere, 
where the cores of the Balmer and He I, 4026, 4471 A lines are 
formed. In this case Teff derived from the He and N ionization 
equilibria should be preferred since the relevant lines form in 
deeper atmospheric layers. 

Since the origin for the discrepancy remains unclear we 
adopted weighted means for the atmospheric parameters and 
had to admit a larger eiTor range for PG 1219+534 than for the 
other programme stars. 

4.2. Abundances 

The resulting abundances are plotted in Fig. ^ and compared 
to PG 1605+072 (paper I). Although KPD 2109+4401 and 
PG 1219+534 do not display C, O, Ne, Mg Al and Si lines, 
stringent upper limits have been derived from the absence of 
the presumedly strongest lines of these species. 

Helium is depleted in all programme stars. Its abundances 
ranges from 1/80 solar (Feige48) to 1/3 solar (PG 1219+534). 

From the interplay of gravitational settling and radiative 
levitation, an anticorrelation between helium abundance and 



gravity might be suspected. However no such correlation ex- 
ists. This has also been found previously from studies of much 
larger samples of non-variable sdB stars (e.g. Schulz et al., 
1991, Fontaine & Chayer, 1997). 

Like helium the metals are deficient with respect to the sun. 
Iron, however, is a notable exception, since it is solar to within 
the error limits in all programme stars. 

Large deficiencies of carbon, oxygen and silicon of up to 
2 dex or more (KPD 2109+4401) have been found, whereas 
nitrogen and sulfur are only mildly deficient. In PG 1219+534 
nitrogen is even almost solar. Large deficiencies of carbon and 
silicon as well as near solar nitrogen abundances have been 
reported for several non-variable sdB stars from high resolution 
UV spectroscopy (e.g. Heber et al. 1984b, Lamontagne et al. 
1985, 1987, see Heber, 1998 for a review). 

It is worthwhile to note that the abundances of C, O, 
Ne, Mg, Al and Si in the high gravity programme stars 
KPD 2109+4401 and PG 1219+534 are considerably lower 
than in the lower gravity stars PG 1605+072 and Feige48 (see 
Fig. ||) which may point to the (selective) action of diffusion, 
i.e. an interplay between gravitational settling and radiative lev- 
itation. It is, however, puzzling that iron is solar irrespective of 
the stellar gravity. 

Recently, Ohl et al. (2000) analysed the far UV spectrum of 
PG 0749+658, a sdB star somewhat cooler than our programme 
stars, obtained with the FUSE satellite and found solar abun- 
dances for several iron group elements. PG 0749+658 shares 
the mild depletion of nitrogen and sulfur with our programme 
stars. KPD 2109+4401 and PG 1219+4401 show a similarly 
strong deficiency of carbon and silicon as PG 0749+658 and 
several other non-variable sdBs (Baschek et al., 1982a,b; Lam- 
ontagne et al., 1985, 1987, Heber et al. 1984b). 

However, until recently only a few diffusion calculations 
were available (e.g. Michaud et al. 1985). Photometric and 
spectroscopic observations of blue horizontal branch stars in 
globular clusters (Grudahl et al., 1999, Moehler et al., 1999, 
Behr et al., 1999a) stimulated calculations of atmospheric dif- 
fusion in atmospheres for such stars (Hui-Bon-Hoa et al., 
2000), which are, however, considerably cooler than our pro- 
gramme stars. The predicted iron abundances are even larger 
than the observed ones. 

For the sdB star PG 0749+658, Ohl et al. (2000) calculated 
equilibrium abundances within the framework of the radiative 
levitation theory and were able to reproduce the observed abun- 
dances of carbon, sulfur and iron, but predicted higher abun- 
dances for nitrogen and silicon and a much lower helium abun- 
dance than observed. Mass loss has frequently been invoked 
to explain the discrepancy between observed and predicted he- 
lium abundances (Fontaine & Chayer, 1997). Unglaub & Bues 
(1998) carried out diffusion calculation for helium, carbon and 
nitrogen for a model of Teff =40 000 K, log g=6.0 including 
mass loss processes. For mass loss rates less than lO^^'^ Mq/yt 
the predicted helium, carbon and nitrogen abundance are quali- 
tatively consistent with the observation. For a quantitative com- 
parison additional calculations at lower effective temperatures 
and gravities as well as for various mass loss rates are required. 
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Fig. 6. Abundances of the programme stars relative to solar values 



Finally, we point out that a solar surface iron abundance is con- 
sistent with the diffusion calculations of Charpinet et al. (1997) 
for the envelopes of pulsating sdB star models. 



4.3. Rotation velocities 

The spectral lines of PG 1605+072 are considerably broad- 
ened, which Heber etal.(1999) attributed to stellar rotation and 
derived vsini = 39km/s, by fitting the strongest metal lines. 
Recently, O' Toole et al. (2000) performed time-series spec- 
troscopy of PG 1605+072 and detected radial velocity vari- 
ations at the same frequencies found from photometry with 
amplitudes up to 14 km/s at H/3. Somewhat smaller veloc- 
ity amplitudes were found for higher Balmer lines. The au- 
thors conclude that these velocity variations arise from the 
pulsations. Since part of the line broadening observed in the 
time averaged spectrum analysed in paper I may be due to ra- 



dial velocity variations, the true projected rotational velocity is 
slightly smaller than vsini = 39km/s. Feige48, PG 1219+534 
and KPD 2109+4401 are all very sharp-lined and we derived 
upper limits of vsini<5-10km/s (see above). This is consis- 
tent with time-series spectroscopy of KPD2109H-4401 (Jeffery 
& Pollaco, 2000) who measured radial velocity variations of 
2 km/s, well below our detection limit. The same holds for any 
rotational broadening that has been inferred from the analy- 
sis of the light curve (Fontaine, 1999, priv. com., Vj-ot sin i= 1.2 
km/s or less). 

The low projected rotational velocities of the pulsating sdB 
stars may be compared to the rotation velocities found for glob- 
ular cluster horizontal branch stars, since sdB stars form the ex- 
tremely hot end of the horizontal branch. Behr et al. (2000a,b) 
analysed Keck HIRES spectra of HB stars in the globular clus- 
ters M 13 and M 15 and found that a star's rotation appears to be 
related to its position along the horizontal branch. Stars hotter 
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Table 7. Radial velocities and equivalent widths of interstel- 
lar Ca II K derived by fitting the observed profile with three 
(PG 1219+534) or four Gaussian profiles. 



Feige 48 


KPD 2 109+4401 


PG 1219+534 


PG 1605+072 


V Wx 


V Wx 


V Wx 


V Wx 


km/s mA 


km/s mA 


km/s mA 


km/s mA 


-55 63 


-42 41 


-57 78 


-55 8 


-38 38 


-13 86 


-28 70 


-39 14 


-17 53 


3 31 


-5 26 


-18 96 


5 12 


11 11 




68 



than about 1 1 000 K were found to be rotating slowly (v sin i < 
lOkm/s) whereas the cooler horizontal branch stars rotate more 
rapidly (vsini « 40km/s for M 13) as previously observed by 
Peterson et al. (1995) for M 13. The low rotational velocities 
of our extreme horizontal branch stars are therefore consistent 
with the slow rotation derived for the hottest stars in M 1 3 and 
M 15. The rotating PG 1605+072 hence appears to be an excep- 
tional case. 

SdB stars are generally believed to evolve directly into 
white dwarfs which are mostly slow rotators as well (Heber 
et al., 1997, Koester et al, 1998). 

5. Appendix: Interstellar Ca II K lines 

Radial velocities and equivalent widths of interstellar Call 
K are derived by fitting the observed profile with three 
(PG 1219+534) or four Gaussian profiles. Results are listed in 
Table 
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